Phase ordering and symmetries of the Potts model 
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We have studied the ordering of the (/-colours Potts model in two dimensions on a square lattice. 
On the basis of our observations we propose that if q is large enough the system is not able to 
break global and local null magnetisation symmetries at zero temperature: when q < 4 the system 
forms domains with a size proportional to the system size while for q > 4 it relaxes towards a non- 
equilibrium phase with energy larger than the ground state energy, in agreement with the previous 
findings of De Oliveira et al. |IlLll2|l . 
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1. INTRODUCTION 

Phase ordering 0, is one of the important topics in non-equilibrium statistical mechanics. For systems with two 
coexisting phases the situation is generally well understood from the analytic and numeric points of view 0, 0, @j . For 
the Potts model [f| with a q > 2-degenerated ground state, the situation is not so clear in general. Also in this case, 
an Allen-Cahn power-law regime and dynamical scaling relations for the structure factor and for the distribution 
of domain sizes 0,H, have been predicted and found numerical ly [Toll. On the other hand, a singular behaviour 
of the Potts model has been recently observed by De Oliveira et. al. when the degeneracy (q) is large enough: 

in the thermodynamic limit the model was shown to relax towards a 'glassy', disordered, phase with a non negligible 
density of defects when it is quenched at zero temperature. 

This paper is devoted to discuss some ideas for understanding this glassy phase. In the next section we recall the 
phase and report new results on the ordering of the Potts model after a quench at low but finite temperature. The 
system equilibrates locally nucleating domains that eventually become of the size of system. In section 3 we study 
the T = case, and we interpret the impossibility to equilibrate in the thermodynamic limit as an impossibility of 
breaking the local zero magnetisation. Section 4 is to conclude and summarise the main results. 

2. DYNAMICS WITH THERMAL FLUCTUATIONS 

De Oliveira et. al. [Til have recently observed an interesting slow relaxation in the dynamics of the Potts model 
after a quench from a disordered state to zero temperature, and the impossibility for the system to achieve the ground 
state in the L —> oo limit, being L the linear size of the system. Let us briefly describe this effect. Given a lattice L, 
in which each site j £ L can take q equivalent values, or colours, the Potts model [f| is defined by the Hamiltonian: 



where {i, j} means i £ L, j is a neighbour of i, and where Cj G {1 • • ■ q} is the colour of site i. When quenching at zero 
temperature a two-dimensional system on a square lattice using Glauber single spin-flip dynamics, systems with q > 4 
obey the Allen-Cahn law 0] e(i) oc i -1 / 2 (e being the energy per site) up to a time i, increasing with L, when they 
get trapped in a blocked configuration, invariant with respect to the single spin-flip dynamics at zero temperature 
|l3l Il4i Il5| |23]|. Inverting the usual order of the t — > oo and L — > oo limits, i.e., assuming that the thermodynamic 
limit is taken before the infinite-time limit, and extrapolating to t — > oo the e(t) data with the Allen-Cahn power-law 
e(t) oc t- 1 ' 2 [13, one finds a positive energy e* = e(t — > oo) > for q > 4 and, in particular the data for e* is 
very well fitted by the expression e* = bL^ 1 / 2 + b'(q — 4) 1 / 2 , b,b' being real constants. This would imply that an 
infinite-size system with q > 4 relaxes after an infinite time towards a phase with stationary observables and positive 
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energy, different from the ground state with zero energy, and this is the reason for which this out-of-equilibrium 
'phase' is called 'glassy phase' in [l^. A noticeable point is that the onset of this glass- like phase occurs in the two- 
dimensional Potts model for q > 4, and this g-range coincides with the one in which the system presents a first-order 
phase transition A justification of this fact is in progress [Tf|. We now give a general characterisation of the 
equilibration of the system at finite temperature, and an interpretation of the non-equilibrium phase in terms of the 
symmetries of the problem in Section 3. 

We have studied the dynamics after a quench at finite temperature, T — 0.1, of the 7 colours Potts model with 
periodic boundary conditions, being T c — ln _1 (l + q 1//2 ) = 0.7730... the critical temperature of the model In Fig. 
1 we present the energy per site (on top) and the magnetisation (centre) with respect to t -1 / 2 of a 2-d square lattice 
system with nearest-neighbours interaction and periodic boundary conditions. The magnetisation, m, is 

™=^(i>2-V9) (2) 

c 

where x c is the fraction of colour c, x c = N c /N, N c being the number of sites in the system with colour c and N = L 2 
is the total number of sites. The magnetisation is zero when x c — Xjq Vc, and one when all spins have the same colour. 
Initially the system is in an uncorrelated configuration at infinite temperature, with e = 2(q — l)/q and m = 0. For 
t < t, systems with q > 4 coarsen with equal proportion of all colours and with the expected power-law dependence of 
the energy on time 6] while, for t > r, finite-size systems equilibrate, breaking the m = symmetry and approaching 
gradually the ground state with m = 1, e = 0. The above description of the problem suggests that, for q > 4 and for 
times larger than r, coalescence effects, not considered in the derivation of the Allen-Cahn law, become relevant in 
the dynamics, and allow the mean domain size, I 0, to grow up to the system size, L. The numerical results reported 
in Fig. 1 suggest that, for T > 0, r is not divergent in L, but constant above a certain L, i.e., that the "nucleation' 
needed to equilibrate the system is a local process no longer dependent on L. In fact, the L = 10 3 and L = 500 
curves of Fig. 1 coincide within a 2 (t,L), the variance of the distribution of energy values corresponding to different 
realisations of the quench (shown in error bars). 

We observe from Fig. 1 that the m = symmetry is broken later in larger systems. As in |17| for the continuous 
Ising model, the magnetisation is supposed to be zero for all times in the absence of finite-size effects, and in fact 
we find 0| a similar scaling relation m(t, L) = m(t/L 2 , 1) for q = 2, 3, and even a slower dependence of m(t, L) on 
L for q > 4. One could ask why r seems to approach a limit value with L —> oo, while the magnetisation at fixed 
times tends to zero in this limit. To answer this question we propose the following argument: systems with q > 4 
equilibrate leaving at t = t the power-law dependence of the energy on time and forming domains that will eventually 
become system-sized, and r seems to be a characteristictime of the model. It is this local nucleation the origin of the 
ordering, rather than the changes on the global colour fraction, which is a finite-size effect |16). In order to supply 
a quantitative support to this argument we define an order parameter, 7, accounti ng fo r the spatial ordering, which 
seems to characterise the dynamical ordering more than the global magnetisation, m [l8| . Let us define 7 as a distance 
from the //(r) = 0, Vr, situation, where //(r) is the magnetisation of a cell centred in r and of size A, independent of 
L, in such a way that X/L — > in the thermodynamic limit. In particular, 

7^i|dr / ,(r) = ^ T ^/dr(^0(r) c 2 -l/ (Z )^-l^U. (3) 

where V is the volume of the general system and ^> c (r) £ [0 : 1] is the proportion of colour c in the cell centered in r. 
The expression at the right of the arrow in (3) is the definition of 7 in the lattice L and cell i is a cell centred in the 
position of site i p4| . For L — > 00, 7 so defined is zero in the completely uncorrelated configuration, and one in the 
ordered configuration, when all the sites have the same colour and m = 1. Since \x > 0, 7 is a distance, functional of 
//(r), and can be used as an alternative order parameter accounting for the spatial ordering of the system. We see in 
Fig. 1 (bottom) that also 7 presents a power-law dependence on time. As expected by the argument exposed above, 
7 seems not to depend on L for large L: curves for L = 500 and 10 3 coincide within their standard deviation limits. 
Moreover, 7 seems to leave the power-law regime and converge to 1 at r. 

3. T — DYNAMICS AND SYMMETRIES OF THE GLASSY PHASE 

At zero temperature, the situation is different: the above defined time, r, needed to leave the Allen-Cahn law is 
divergent with L, as described in [Tl|. and an infinite-size system is always supposed to follow the Allen-Cahn power- 
law. In larger and larger systems, it is less and less likely to find a blocked configuration or a path in phase space to 
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the ground state. A measure of this fact can be seen in the value of the energy variance, a 2 (t,L), which decreases 
with L (see inset of Fig. 3). This means that in a large system it is less likely to find a deviation from the power-law 
regime, in which blocked dynamics ^| or breaking of the m = symmetry [T^ are not present. This explains at a 
qualitative level the fact that in large systems the magnetisation is broken later (and hence the system equilibrates 
later), but it does not address the fact that the energy of an infinite size system with q > 4 converges to a positive 
value for t — > oo. In other words, it does not explain why the term e* of the generalised Allen-Cahn law e = at^ 1 / 2 + e* 
is different from zero for q > 4. One could ask why systems with q > 4 cannot converge to a zero energy phase in the 
limit L — * oo, even respecting the m = symmetry, as the q = 2, 3 cases do. We discuss this point in the following. 

Numerical simulations show that in the latest stage of the coarsening at zero temperature, systems with q = 2, 
3 and with large enough L (to respect the m = symmetry) present final configurations formed by domains of the 
size of system (Fig. 2, left) and i(t), the mean domain size, is proportional to L when t — > oo. The energy, or the 
perimeter of the interface, is proportional to L. and thus, the energy per site is zero in the thermodynamic limit. In 
fact we have e*(q) = for q = 2, 3, as said before. On the other hand, for q > 4, we observe that the system presents 
not only m — 0, but also a local symmetry that is a local equal fraction of all colours. We define this local symmetry 
by introducing a certain scale, A, with L > A > t, such that the magnetisation is also zero in every cell of the system 
of size greater than A and being A such that X/L — > for L — > oo. We propose, and we argue below, that the zero 
temperature dynamics cannot form domains of the size of the system for q > 4, hence the perimeter of the interface 
(as well as the final number of domains) grows with L 2 and the energy per site of the glassy phase is not zero in the 
thermodynamic limit, while, in this limit, systems with q = 2, 3 do not break the m = symmetry when forming 
domains of the size of the system but they can break the local m = symmetry. A numerical confirmation of this 
argument is that 7 converges to 7* = 1 in the L, t — ► 00 limits for q = 2, 3, but 7* < 1 for q = 7: the fraction of the 
volume V for which fi ^ 1 is not negligible in the thermodynamic limit. In Fig. 3 we report the energy, magnetisation 
and order parameter of local ordering as a function of t~ x l 2 for the T = dynamics in the q — 7, 3 and 2 cases. It 
seems that 7 presents the same power-law dependence on time as the energy in the q = 7 case. Inverting the order 
of the limits and extrapolating to t — + 00 with a linear fit, as done for the energy in 111, we have 7* = 0.807 ± 0.005, 
while for q = 2, 3, 7 converges to 1. 

This local zero magnetisation for q > 4 (see Fig 2, right) characterises the glassy phase also in the sense that it 
predicts a q 1 / 2 dependence of e*, which coincides in the q — ► 00 limit with the result e* ~ (q — 4) 1 / 2 of 01 . Let us 
show this fact in the approximation of domains of identical size £ |19| • In order to satisfy the local-m = condition 
there must be nq domains inside each cell of size A, n being an integer. In d dimensions it is e ~ nr> l d ~ l jL d = , 
where no = L d £~ d is the number of domains, and setting \ d ji d = nq we have e ~ \~ 1 q 1 / d . When q approaches 4, 
coalescence effects not considered in this argument become important, to which we attribute the slower dependence 
on q of the energy, e* ~ (q — 4) 1 / 2 . 

As we have seen there are evidences that, in absence of thermal fluctuations, the Metropolis dynamics cannot break 
the local symmetry of colours when q > 4. This fact and the q > 4 limit can be explained (not yet in a rigorous 
fashion) as a consequence of an assumption on the glassy phase symmetries. To introduce it let us define the q 
'coarse-grained' fields <^> c (r) G [0, 1], in analogy with as the fraction of colour c in a cell centred in r and of size 
A. Clearly J^t fic = 1 and the fields <f> c describe the configuration in a A-dependent way. Given the fact that the 
Hamiltonian is invariant with respect to S q , the group of permutations of the q colours, we assume that the m — 
symmetric configuration corresponding to L — > 00, t — > 00 is maximally ordered in such a way that any permutation 
of the colours, p € S q , is equivalent to a p-dependent global spatial transformation, Up. i.e.: 

<t> p(c) (r) = MU P r) VpeS q (4) 

If we assume the operator U to be linear we have that U : S q — > jSf (R 2 ) is a representation of S q in the vector 
space of linear operators in the plane, Jz? (M 2 ). Since there is no representation of S q in ££ (R 2 ) for q > 4 different 
from the trivial representation (U p — I, the identity operator on ffi 2 , Vp), from the above assumption it follows the 
local-zero magnetisation of the glassy phase for q > 4: setting U p = I in (4) we obtain 4> c i{r) = 4> c (r), that with the 
normalisation condition gives 4> c ( r ) = 1/? Vr, c, i.e., the local m — symmetry in the cell of size A. The assumption 
(4) can also 'justify' the spatial symmetries of configurations with q < 4 in the m = regime: for q = 2 the group S2 
admits a nontrivial representation in the plane, {I, a}, where I is the identity and a is the rotation of ir radians in 
the polar angle and, in fact, configurations with q = 2 in the glassy phase seem to present this symmetry (see Fig. 2, 
left), and a similar 27r/3 rotational symmetry seem to exist for q — 3 configurations, even if it is often hindered by the 
two-preferred directions of the square lattice interaction 0] and by the formation of structures locally stable under 
the T = dynamics fTL l2Qj - Moreover, in a general d-dimensional system, there exist no nontrivial representations of 
S q in M d for q > d + 2 [22j . This proposition seems to be true for all q > 4, and we are workingon the general proof. 
If verified, this consequence of the assumption (4) would coincide with the result by Lifshitz [lj, who argued that a 
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d-dimensional system quenched below its critical point does not necessarily equilibrate into an ordered state, in the 
presence of a ground state that is degenerated more than d + 1 times. 



4. CONCLUSIONS AND FURTHER RESEARCH 



We have studied the ordering dynamics of the 2-d Potts model in a square lattice with q = 2, 3 and 7 by use of 
single-spin flip dynamic Monte-Carlo simulations. At positive temperature, T — 0.1, and for q = 7 the system relaxes 
leaving the Allen-Cahn power law at a given time, which is independent on L for large L. At zero temperature, in 
the L — > oo limit, systems with q > 4 are not able to nucleate breaking the local symmetry, and hence they converge 
to a phase with nonzero energy. We give a quantitative support for this argument in the q = 7 case with the help of 
an order parameter of local ordering that we define. The q > 4 limit is presented as a consequence of a hypothesis 
on the symmetries of the glassy phase. For both temperatures the magnetisation at equal times decreases with 
increasing sizes and it is argued to be zero in the thermodynamic limit at any time \l7\ - A generalisation of this 
study for different q values and for larger L is being performed [l6l ] together with an investigation of the relation 
between the ordering properties and the phase space structure. 
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FIG. 1: Energy per site (e), magnetisation (m) and order parameter of local ordering (7) with respect to t™ 1 / 2 for the Potts model 
with 9 = 7 after a quench to temperature T = 0.1. Results are averaged over 200, 60 and 30 realisations (for L = 200, 500 and 10 3 , 
respectively). Time is in MCS units. Error bars in the i/-axes are the variances, <r 2 , of e and 7, corresponding to the average over different 
realisations of the quench (shown every 2000 MCS only). In larger systems the m = symmetry is broken later Il7l . c 2 decreases with 
the size of the system at equal times. The blue line is a fit to the L = 10 3 system energy in the range t -1 / 2 £ [0.02, 0.1], and the vertical 
dotted line is an estimation of the nucleation point, corresponding to the time at which e differs from the fit more than cr 2 /2. At that 
time, t, the system nucleates leaving the Allen-Cahn power law, and 7 increases towards 1. 




FIG. 2: Example of two configurations reached by the 2 and 7-colours Potts in a L = 500 lattice, after 6.25 10 4 MCS of a single realisation 
of a quench at zero temperature. The magnetisation is m = 1.9 10 -4 and 9.6 10 -4 , respectively. Spins with a given colour are shown in 
blue, while spins with the remaining q — 1 colours are shown in white. In the latest stage of the coarsening the interfaces between domains 
of the q = 2 case become straight lines. 
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FIG. 3: Energy per site (e), magnetisation (m) and order parameter of local ordering (7) with respect to t~ 1 ' 2 for the Potts model after 
a quench at zero temperature. The q = 7, 3 and 2 cases are presented, and three different sizes for the q = 7 case. Magnitudes are an 
average over 240, 60 and 40 realisations of the quench in systems with L = 200, 500 and 10 3 , respectively . Fo r the cases q = 2 and 3 the 
energy converges to zero following the Allen-Cahn law, while in the q = 7 case it converges to e* = 0.1192 The global magnetisation 

decreases with q at equal times and sizes, for all sizes and times studied. It also decreases with L for all times |17|. and it is supposed to 
be zero in the thermodynamic limit. Even with this constraint, systems with q < 4 can form system-sized domains with zero energy, and 
in fact they order spatially, and -y(t — > 00) — > 1, while the zero-temperature dynamics cannot break the local zero magnetisation symmetry 
to form domains of the size of the system when q > 4, as illustrated by the fact that -y(t — » 00) — » 7* < 1 for q = 7, and this is the 
reason for which e* > 0. The extrapolated line is a fit to the L = 10 3 data for 0.005 < t^ 1 / 2 < 0.06 MCS -1 / 2 . The fit coincides with the 
data within its variance limits. The equation of the fit is 7f it = 7* + at -1 / 2 with 7* = 0.8072 and a = —3.1756. On inset we present the 
variance <r 2 of the energy, which is a decreasing quantity with L. for all times and for the three sizes studied. 



